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Numerical  and  experimental  investigation  of  the 
turbulent  flow  in  a ribbed  serpentine  passage 

By  Giemluca  laccarino,  Georgi  Kalitzinf  and  Christopher  J.  Elkins| 


1.  Introduction  and  motivations 

Modern  gas  turbine  engines  operate  at  high  combustor  outlet  temperatures  to  achieve 
higher  thermal  efficiency  and  thrust.  Turbine  blades  are  exposed  to  these  high-temperature 
gases  and  undergo  severe  thermal  stress  and  fatigue.  The  design  of  highly  efficient  cool- 
ing systems  for  turbine  blades  has  an  enormous  potential  impact  on  engine  development. 
Cooling  devices  are  based  on  a secondary  flow  system  built  into  each  blade,  as  illus- 
trated in  Fig.  1.  The  secondary  flow  passages  are  extremely  complicated  consisting  of 
one  or  multiple  legs  with  turbulators  (rib-roughened  serpentines),  holes  connecting  the 
secondary  path  to  the  external  surface  of  the  blade  (film  cooling),  tube  bundles,  slots, 
etc.  The  geometrical  complexity  of  these  passages,  however,  is  extremely  challenging  for 
the  use  of  advanced  simulation  tools  based  on  state-of-the-art  three-dimensional  CFD 
solvers.  The  generation  of  a computational  grid  (even  using  unstructured  mesh  tech- 
nology) requires  a considerable  amount  of  time.  As  a result,  the  analysis  of  the  cooling 
performance  of  the  system  is  largely  based  on  isolated  sub-component  simulations,  sim- 
plified one-dimensional  models  and  experimental  correlations. 

A large  amount  of  research  has  been  devoted  to  the  analysis  of  the  flow  characteristics 
(in  particular  tmbulence  statistics)  and  heat  transfer  in  channels  with  ribs.  The  effects 
of  rib  shape,  pitch,  height  and  inclination  with  respect  to  the  incoming  flow  have  been 
quantified  from  an  experimental  point  of  view.  Numerical  predictions  of  these  flows  are 
complicated  because  of  the  increased  turbulence  intensity  near  the  ribs  and  its  effect  on 
the  flow.  In  Reynolds-Averaged  Navier-Stokes  (RANS)  approaches  the  turbulence  inten- 
sity is  strongly  affected  by  the  choice  of  the  turbulence  closure.  RANS  eddy  viscosity 
models  based  on  the  linear  Boussinesq  relationship  between  Reynolds  stresses  and  ve- 
locity gradients  are  the  de  facto  standard  in  industrial  numerical  simulations.  Previous 
investigations  (Ooi  et.  al.  2002a;  lacovides  & Raisee  1999)  have  illustrated  the  predic- 
tive capabilities  of  such  models  in  turbomachinery  cooling  systems  components:  smooth 
channels,  channels  with  ribs,  180  degrees  turns  (U-bends),  etc.  In  particular,  it  has  been 
shown  that  the  flow  is  strongly  three-dimensional  and  the  presence  of  oblique  ribs  gen- 
erates strong  coherent  vortices  in  the  domain  (Ooi  et  al.  2002b).  It  must  be  noted  that 
most  of  these  simulations  are  performed  under  the  assumption  that  the  flow  is  periodic  in 
the  streamwise  direction  and,  therefore,  only  a small  section  of  the  passage  is  effectively 
simulated. 

Several  experimental  datasets  are  available  to  evaluate  the  accuracy  of  these  simu- 
lations (Han  et.  al.  1985;  Rau  et.  al.  1998);  typically  the  measurements  are  taken  in 
devices  that  include  many  ribs  in  between  duct  U-bends  and  the  data  axe  collected  at 
different  locations  to  verify  that  periodicity  conditions  are  indeed  applicable.  Although 
these  comparisons  between  experiments  and  computations  are  invaluable  in  improving 
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Figure  1.  Sketch  of  the  Cooling  System  of  a Turbine  Blade  (Cumpsty  1997). 

the  understanding  of  the  physics  and  in  providing  accurate  model  validation,  they  do  not 
represent  the  wealth  of  complexity  and  interactions  present  in  more  realistic  configura- 
tions, like  a multiple-leg  serpentine. 

Unfortunately  the  computational  cost  of  such  simulations  becomes  substantial.  For  a 
single  three-dimensional  rib-roughened  passage  it  might  be  necessary  to  generate  a grid 
with  half  million  cells  (Ooi  et.  al.  2002b)  to  completely  resolve  the  flow  structures.  In  a 
serpentine  with  thirty  ribs  and  two  180  degree  turns,  the  mesh  size  should  be  in  the  order 
of  30  million  cells.  In  addition,  it  is  necessary  to  obtain  extensive  experimental  datasets 
to  quantify  clearly  the  accuracy  of  the  simulations. 

It  is  clear  that  a novel  numerical  and  experimental  technique  for  investigating  flows  in 
complex  geometries  would  be  extremely  beneficial  to  the  development  and  improvement 
of  the  predictive  tools  used  in  the  design. 

In  this  paper,  the  turbulent  flow  in  a serpentine  with  oblique  ribs  is  investigated  ex- 
perimentally and  by  numerical  simulations.  The  measurements  are  carried  out  by  using 
Magnetic  Resonance  Velocimetry  (MRV)  (Elkins  et.  al.  2003)  and  the  simulations  using 
the  Immersed  Boundary  (IB)  technique  (laccarino  & Verzicco  2003;  Kalitzin  & laccarino 
2002).  A brief  description  of  these  two  approaches  is  reported  in  following  sections.  The 
results  are  reported  in  terms  of  velocity  distributions  in  various  planes  in  the  serpen- 
tine; differences  between  measurements  and  simulations  are  presented  qualitatively  and 
quantitatively.  The  study  of  the  discrepancy  allows  us  to  identify  areas  of  needed  im- 
provements in  the  turbulence  modeling. 


2.  Ribbed  serpentine 

The  serpentine  used  in  this  study  is  represented  in  Fig.  2;  the  model  was  drawn  in 
AutoCAD  and  fabricated  using  a 3D  System  250-50  SLA  machine.  The  serpentine  has  a 
square  cross  section  of  height  H = 20mm  and  ten  staggered  oblique  ribs  on  the  top  and 
bottom  walls.  The  rib  height  is  Q.IH  and  the  pitch  (distance  between  two  successive  ribs) 
is  O.Giif;  the  rib  angle  is  45  degrees.  A fully  developed  pipe  flow  enters  the  first  leg  of  the 
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Section  A-A  Section  B-B 


Figure  2.  Sketch  of  the  cooling  passage  investigated  and  the  region  where  the  experiments 
and  computations  have  been  compared. 


serpentine  through  a converging  section  yielding  a uniform  velocity  profile.  Two  Reynolds 
numbers  have  been  investigated,  8,000  and  50,000  based  on  the  serpentine  height  and 
the  bulk  velocity;  in  the  present  paper  only  the  results  at  the  low  Reynolds  number  are 
reported. 


3.  Experimental  technique 

Magnetic  Resonance  Velocimetry,  MRV,  has  been  used  to  measure  the  mean  velocity 
components  in  the  entire  serpentine.  A massive  amount  of  data,  corresponding  to  about 
2 million  locations  with  1mm  resolution  (0.05H)  has  been  collected.  The  measurement 
technique  is  described  in  detail  in  Elkins  et.  al.  (2003)  where  a simple  pipe  flow  and  a 
serpentine  similar  to  the  one  studied  here  (Fig.  2)  are  investigated. 

The  MRV  technique  is  a non-invasive  experimental  method  for  measuring  mean  veloc- 
ities and  is  implemented  in  medical  magnetic  resonance  imaging  systems.  MRV  has  been 
used  for  many  low-Reynolds  number  investigations  (Fukushima  1999),  such  as  physio- 
logic, multiphase,  and  porous  media  flows.  It  has  also  been  used  in  a few  simple  turbulent 
flows  such  as  pipe  and  jet  flow.  Our  current  focus  is  in  high  Reynolds  number  turbulent 
flows  in  complex  internal  geometries  with  the  objective  of  measuring  the  time-averaged 
three-component  velocity  field  in  order  to  provide  physical  insights  and  an  extensive 
database  for  CFD  validation.  It  should  be  mentioned  that  the  MRV  technique  is  com- 
patible with  the  rapid  prototyping  processes  of  stereolithography  and  fused  deposition 
marbining.  Utilizing  current  RP  technology  highly  complex  flow  models  can  be  designed 
and  manufactured  in  only  a few  days,  and  the  MRV  measurements,  typically  producing 
over  1 million  velocity  vectors,  can  be  completed  in  30  minutes. 
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Figure  3.  Resolution  used  in  the  coarse  grid  simulations  (left)  and  experiments  (right)  in  the 

region  5 (see  Fig.  2). 


4.  Numerical  technique 

The  numerical  simulations  presented  in  this  work  are  based  on  the  Immersed  Boundary 
(IB)  technique.  The  method  is  presented  in  detail  in  laccarino  & Verzicco  (2003)  and 
Kalitzin  & laccarino  (2002);  only  a brief  description  is  reported  herein. 

In  the  IB  method  Cartesian  non-uniform  meshes  are  used  and  forcing  terms  are  added 
to  the  governing  equations  to  enforce  boundary  conditions  (no  slip  velocity,  for  example) 
on  surfaces  not  aligned  with  the  grid  lines.  The  computational  code  used  is  based  on  the 
steady-state  incompressible  RANS  equations  closed  with  either  a one-equation  (Spalart 
k Allmaras  1994)  or  a two-equation  model  (Kalitzin  & laccarino  2002).  A second-order 
implicit  discretization  is  used  within  a SIMPLE  pressure-velocity  coupling  algorithm;  the 
turbulence  equations  are  solved  in  a segregated  manner.  Multigrid  is  used  to  accelerate 
the  convergence.  The  IB  forcing  is  based  on  a tri-linear  reconstruction  (laccarino  & 
Verzicco  2003)  for  the  velocity  components  and  the  turbulent  scalars. 

The  geometry  definition  is  based  on  the  STL  file  format  and,  therefore,  the  same 
AutoCAD  model  used  to  build  the  experimental  configuration  has  been  used  to  obtain  a 
simplified  CFD  computational  model.  The  computational  grid  and  the  IB  interpolation 
stencil  is  generated  automatically  and  the  pre-processing  time  to  start  the  computation 
is  negligible. 


5.  Results  and  discussion 

The  geometrical  configuration  is  shown  in  Fig.  2.  The  CFD  domain  includes  only  the 
first  three  legs  of  the  serpentine.  The  simulations  have  been  performed  using  computa- 
tional grids  ranging  from  half  million  to  forty  millions  elements.  An  example  of  the  coarse 
grid  resolution  of  the  CFD  mesh  is  shown  in  Fig.  3.  The  underlying  grid  for  the  MRV 
experiments  is  shown  in  the  same  figure. 

The  fine  grid  simulations  are  compared  to  the  experimental  data  in  the  areas  indicated 
in  the  sketch  of  Fig.  2.  In  Fig.  4 and  5 two  horizontal  planes  are  reported  midway 
through  the  ribs  and  halfway  through  the  passage,  respectively.  In  Fig.  7 and  8 two 
vertical  sections  are  reported  with  the  streamwise  direction  pointing  away  from  the  page 
and  toward  it,  respectively.  Finally,  in  Fig.  6 a vertical  longitudinal  plane  showing  the 
streamwise  flow  (from  left  to  right)  is  reported. 


0.0  0.1  0.4  OJ  0.7  0.9  l.l  1.3  1.4  l.< 

velocity  magnitude 


Figure  5.  Comparison  between  simulations  (top)  and  experiments  (bottom)  in  region  2 (see 
Fig.  2):.  horizontal  plane  corresponding  to  a cut  at  the  mid  height  of  the  passage. 


The  general  qualitative  agreement  between  experiments  and  simulations  is  remarkable. 
In  particular,  the  secondary  flows  in  between  the  ribs  described  in  Ooi  et.  al.  (2002b) 
are  well  represented.  The  flow  in  between  the  U-bends  shows  a repeating  patterns  after 
about  three  or  four  ribs  demonstrating  the  importance  of  the  complete  simulation  as 
opposed  to  a periodic-rib  channel  study. 

The  simulations  capture  the  most  dominant  effects  illustrated  in  the  experiments: 

(o)  A strong  three-dimensionality,  as  indicated  by  the  mid-height  plane.  High  velocity 
regions  oblique  to  the  direction  of  the  flow  (and  orthogonal  to  the  rib  orientations) 
indicate  the  presence  of  strong  streamwise  vortices; 

{h)  Recirculation  bubbles  in  the  downstream  region  of  the  180  degree  bends;  in  partic- 
ular these  features  are  in  remarkable  good  agreement  (in  terms  of  length  and  maximum 
reverse  velocity)  with  the  measurements; 

(c)  A massive  separation  behind  each  rib.  Fig.  7.  Again  the  length  and  the  velocity 
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x-velocity  component  y- velocity  component  z-velocity  component 


Figure  6.  Comparison  between  simulations  (top)  and  experiments  (bottom)  in  region  5 (see 
Fig.  2):  vertical  plane  corresponding  to  a longitudinal  section  before  the  U-bend. 

magnitude  are  in  good  agreement,  even  if  the  resolution  in  the  experiments  is  limited 
(cfn.  Fig.  3). 

(d)  High  speed  regions  atop  of  each  rib  with  high  shear  and  turbulence  (not  reported). 

In  addition,  they  clearly  indicate  that  the  flow  in  the  straight  sections  is  accurately 
represented  and  the  largest  discrepancies  are  around  the  U-bend  regions.  This  is  a well 
known  limitation  of  the  (linear)  eddy  viscosity  models  used  in  the  RANS  simulations  and 
it  does  not  appear  to  be  linked  to  the  Cartesian  nature  of  the  flow  solver. 

The  wealth  of  information  provided  by  the  MRV  technique  allows  a clear  identification 
of  the  accuracy  of  the  simulation.  In  particular,  it  can  be  used  to  demonstrate  the  limi- 
tation imposed  by  the  turbulence  model.  The  three  components  of  velocity  are  reported 
for  various  flow  regions  in  Fig.  6,  7 and  8.  Region  5 and  3 are  located  before  the  U-bend. 
The  flow  is  here  fully  developed  and  the  streamwise  component  in  the  cross  section  shows 
a typical  inverted  C shape,  with  larger  velocities  near  the  left  wall.  This  combined  with 
the  plot  of  the  cross-flow  component  (y-velocity,  cfh.  Fig.  2)  shows  the  presence  of  two 
fairly  symmetric  (with  respect  to  the  channel  mid-height)  longitudinal  vortices.  On  the 
other  hand.  Fig.  8 shows  the  same  cross-section  after  the  U-bend,  and  in  this  case  the 
vortical  structures  are  not  yet  formed  and  the  flow  is  substantially  affected  by  the  bend. 
In  this  section  the  difference  between  experiments  and  simulations  appears  to  be  larger, 
suggesting  a limitation  in  modeling  turbulence  in  the  bend  region. 

Finally,  the  error  in  the  computed  velocities  is  reported  in  Fig.  9 for  a vertical  section 
before  the  U-bend.  The  CFD  results  are  interpolated  on  the  experimental  locations  and 
the  differences  in  the  three  velocities  components  (normalized  by  the  experimental  values) 
are  plotted.  The  largest  error  is  about  5%  and,  most  importantly,  has  a unorganized 
pattern,  showing  that  the  same  flow  features  (streamwise  vortices  and  secondary  rib 
structures)  are  present. 
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Figure  7.  Comparison  between  simulations  (top)  and  experiments  (bottom)  in  region  3 (see 
Fig.  2):  vertical  plane  corresponding  to  a cross-section  before  the  U-bend. 


Results  obtained  on  coarse  grids  show  the  same  main  features  and  agree  fairly  well 
with  the  experiments. 


6.  Conclusions  and  future  work 

A numerical  and  experimental  investigation  of  the  turbulent  flow  inside  a multiple-leg 
serpentine  with  oblique  ribs  is  presented.  Both  the  simulations  and  the  measurements  are 
based  on  the  use  of  novel  techniques:  the  immersed  boundary  approach  for  the  numerical 
predictions  and  the  magnetic  resonance  velocimetry  for  the  experiments.  Both  technique 
are  very  fast.  Once  the  experimental  stand  is  setup,  large  amount  of  data  can  be  measured 
in  very  short  time  for  a CAD  model  that  fits  into  the  measuring  section.  The  same  CAD 
model  is  used  to  set  up  the  computations  without  extra  cost  in  time. 

Extremely  detailed  comparisons  axe  reported  due  to  the  high  resolution  of  the  two 
approaches  considered:  up  to  40  million  grid  cells  are  used  in  the  simulations  and  about 
2 million  point  measurements  are  available.  The  results  show  an  overall  satisfactory 
agreement  and  the  turbulence  model  used  (two-equation  model)  seems  to  perform  quite 
well  for  the  straight  parts  of  the  passage  capturing  the  three-dimensional  flow  structures 
accurately.  The  major  discrepancies  are  present  in  the  U-bend  and  in  the  immediate 
downstream  areas  where  the  linear  eddy  viscosity  approaches  used  cannot  account  for 
the  strong  turbulence  anisotropy  (lacovides  & Raisee  1999). 

Additional  simulations  and  experiments  have  been  carried  out  at  a larger  Reynolds 
number  and  the  results  show  again  an  encouraging  agreement.  Future  work  will  be  de- 
voted to  perform  similar  comparisons  between  these  two  techniques  in  different  complex 
configurations;  in  particular,  in  the  cardiovascular  system. 
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Figure  8.  Comparison  between  simulations  (top)  and  experiments  (bottom)  in  region  4 (see 
Fig.  2):  vertical  plane  corresponding  to  a cross-section  after  the  U-bend. 
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Figure  9.  Differences  between  computations  and  experiments  in  region  3 (see  Fig.  2):  vertical 
plane  corresponding  to  a cross-section  before  the  U-bend. 
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